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Abstract. We provide a renormalization procedure for ^-derivable approxima- 
tions in theories coupling different types of fields. We illustrate our approach 
on a scalar (p^ theory coupled to fermions via a Yukawa-like interaction. The 
non-perturbative renormalization amounts to fixing the scalar coupling via a 
set of nested Bethe-Salpeter equations coupling fermions to scalars. 
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1. Introduction 

The two-particle-irreducible (2PI) effective action has recently regained interest 
in many fields of research such as the QCD phase diagram [ or the dynamics 
of quantum fields out-of-equilibrium [|3j. There are however two major difficulties 
related to how fundamental properties such as UV finiteness or gauge invariance 
manifest themselves within a given truncation of the 2PI effective action. Recent 
progresses have been made concerning renormalization in the case of scalar theories 
[0]. Aside from the question of gauge invariance, an important step towards the 
understanding of renormalization in the case of gauge theories consists in considering 
a theory coupling multiple fields. 

As an illustration, we consider here a massless fermionic field ip at finite tem- 
perature, coupled to a massless self-interacting scalar field ip {Xip'^ theory) via a 
Yukawa-like interaction gipipip. We concentrate on the symmetric phase for which 
it is enough to consider the 2PI effective action as a functional of the full propagators 
S and D: 



[S,D] = - 



^S,D], (1) 
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where $ is the sum of 0-leg skeletons (2PI diagrams). When evaluated at its station- 
ary point, the 2PI effective action provides a powerful way to compute the pressure 
of the system. Further thermodynamic quantities are obtained by taking deriva- 
tives. The propagators are given by the stationary condition on FapJS', D] which is 
suitably translated into a set of equations of motion for the self-energies S and 11. 
As an illustration we consider here the one-loop approximation; 

= -g'i D{Q)S{Q + P) -IfSZ^ , 

JQ 

Il{K) = ^{\ + SX)^^D{Q)+ghr^^S{Q)S{Q + K)+K^SZ^. (2) 

The sum- integrals in the above equations are UV divergent. Renormalization con- 
sists in building temperature independent counterterms SZ^, dZ^ and SX which 
absorb all the divergences. In dimensional regularization there is no need for mass 
counterterms. At one-loop level, the Yukawa coupling is not renormalized and we 
thus take 6g = 0. A general discussion to higher loop orders may be found in [l^j. 

2. Renormalization 

One starts focusing on the zero temperature solutions to the equations of motion 
that we denote by 5x^0 and £'t=o- The related self-energies are denoted by St=o 
and IIt^o and determine the field strength counterterms SZ^ and SZ^ through the 
renormalization conditions dI]T=o/dl^\p^=~fj,^ = and dIlT=o/dK^\K2^_^2 = 0. 

The next step consists in removing coupling singularities from the temperature 
dependent contributions St = S — St=o and IIt = 11 — IIt^o- It is easy to 
check that the equation for St does not contain any UV divergence. In contrast 
the equation for IIt contains logarithmic divergences. A diagrammatic analysis [ 
1^ reveals that these are exactly the same than those encoded in the four-point 
function with four scalar legs (T^^). Renormalization of these divergences thus 
amounts to imposing the renormalization condition r^^(P*,i4r,) = A. To proceed, 
we thus need to know how to build T^^, show that it can be renormalized and check 
that its renormalization removes the coupling divergences in IIt- 

The function is built in two steps from a set of nested Bethe-Salpeter 
equations. These equations involve 2PI kernels which are directly related to deriva- 
tives of the functional $ namely A^^ = -~S^^/SS6S, A^^ = ~2S^^/SSSD, A^^ = 
~2S'^^/dDSS*' and A^^ = AS^^/SD6D. In the first step, one builds up the four- 
point function with four fermionic legs T^^ from A^^: 

r^4P,K)^A^^iP,K)- [ A^^{P,Q)M{Q)T^^iQ,K). (3) 

Jq 

is then combined with the 2PI kernels in order to generate a new kernel: 

A^^iP,K) = A^^{P,K)+ [ A^4P,Q)MiQ)A^^{Q,K) 

Jq 
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A^^{P,Q)M{Q)T^4Q, R)M{R)A^^{R, K) . (4) 



The explicit distribution of fermionic indices in Eqs. Q and 10} can be found in [ 
|Sj where we also define M as M(^a0},{f5){Q) = Saj{Q)Ssi3iQ)- In the second step, a 
second Bethe-Salpeter equation builds up from A^p^p: 

V^p{P,K)=K^^[P,K)~ [ K^^[P,Q)D\Q)V^^{Q,K). (5) 

JQ 

At one loop level, it is simple to check that equation Q is finite. In contrast, 
Eq. Q contains logarithmic divergences. However these are overall divergences 
from which it follows that differences such that h.^ip{P, Q) — h.^ip{K, Q) are finite. 
Furthermore, by construction, the kernel A^^ is 2PI with respect to scalar lines 
implying that A^^{P, Q) — A^^{K, Q) ^ 1/Q at large Q and fixed P and K It 
is then possible to rewrite the equation for r^i^ in a UV finite form: 

T^^iP,K)-T^^{P,,K,) - A^^iP,K)-A^^iP,,K,) (6) 

1 



{Ayy(P, Q) - Ap^(P.,Q)}d^{Q)T^^{Q, K) 
T^^{P,,Q)D^{Q)[a^p{Q, K) - A^^iQ, if*)} . 



Using that T^^{Q,K) - logQ, D{Q) - and A^^{Q, K) ~ A^^{Q, K,) - 1/Q 

at large Q, one checks that the integrals in ijHJl are UV finite. Thus T^^ is finite, as 
announced. 

Finally, in order to show that renormalization of F^ip simultaneously removes 
the coupling divergences in the equation for IIt, it is convenient to express the 
latter in terms of A^^. FoUowings the steps presented in [ 5 , one obtains: 

IIt{K) = ^ f A^p{K,Q)SD{Q) + ^ f A^^iK,Q)a4Q)+nr{K), (7) 

J Q J Q 

where we have performed the Matsubara sum in order to separate implicit and 
explicit thermal dependences. The first integral in (|7|) involves 6D{Q) = D{Q) — 
-Dt=o(Q) = -nTiQ)D^{g) + Dr{Q) with Dr{Q) ~ l/Q^ at large Q. The sec- 
ond integral involves <Tip{Q) — e(go)?^(ko|)P¥>('?o, 9)j a particular combination of the 
sign function e(go), the scalar thermal factor n(|(7o|) and the scalar spectral density 
Pip{qo, q). Q designates integration along the real axis in contrast to Q which desig- 
nates integration along the imaginary axis. Finally Ilr{K) is a finite function which 
decreases as at large K. Using Eqs. © and Q one shows that 

nT{K)^Ilr{K)+^ J_T^^{K,Q)apiQ) + ^ J T^^iK,Q){Dr{Q)^D\Q)Ilr{Q)} . 

(8) 

Using the finiteness of Tip^p and the asymptotic properties of Tip^p, cr^p, Dr and 11^, 
it is simple to check that this last equation is finite, as expected. 
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3. Conclusions 

We have shown on a particular example how to implement renormalization of <i>- 
derivable approximations in theories with multiple fields (see also [ 5 ). This work 
is particularly important for gauge theories where one needs to disentangle the UV 
divergences related to gauge, matter and ghost fields. 
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